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v  ABSTRACT  Cg-Cxe’V^Y^  t  ti<^\i<kOt) J 

^^The  combination  of  ab  initio  calculation  of  the  electronic  wavefunctions  with  a 
wave  packet  calculation  of  the  nuclear  motion  is  used,  within^  the  Born-Oppenheimer 
approximation,  to  compute  the  vibrational  and  electronic  absorption  of  a  polyatomic 
molecule.  A  particular  virtue  of  this  approach  is  that  high  as  Well  as  low  temperature 
spectra  are  both  calculable.  This  method  is  applied  to  C^H,  for  which  the  Complete 
Active  Space  Self  Consistent  Field  1CASSCF)  method  is  used  to^jietermine  full  Born 
Oppenheimer  potential  surfaces.  Using_lhe  assumption_Lhat  the  A(2n)<-  X(2£  +  )  absorp¬ 
tion  can  be  written  as  the  sum  of  the  A(*a)<-  X  and  /(^A")*-  X  absorptions,  the  spectra 
are  determined  to  6QcriT^ resolution  at  a  temperature  of  3000K.  As  a  result  of  the  large 
thermal  bending  amplitude  at  3000K,  the  calculated  spectra  are~-biqad  and  have  little 
resolved  structure.  Two  bands  are  resolvable,  one  is  due  tp.  the  A(*a  )f-X  .absorption 
jind  is  centered  at  SSOQcrn21^  while  the  other  is  doe  to-^he  A(‘A')<-  X  absorption  and  is 
centered  at  950Qcmu*T  The  dramatic  blue  shift  of  th/e  A(*A )—  X  band  results  from  the 
combination  of  the  large  X  state  thermal  bending  amplitude  and  the  high  bending  fre- 
I  auency  of  the  Af2^")  state.  We  also  determine  the  X  state  pure  vibrational  absorption 
^  spectrum  and  show  it_Lo  be  of  much  lower  intehsitv  than  the  mire  electronic  snertrnm  merr.. 
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.spectrum  and  show  it  xti  be  of  much  lower  intepsity  than  the  pure  electronic  spectrum-vst. 
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I.  introduction  4W  r  ft  t^-Vj 

Small  polyatomic  molecules  are  often  found  in  high  temperature  environments  such 
as  flames,  stars,  and  shock  fronts.  From  a  knowledge  of  their  cross-sections,  we  can 
learn  much  about  the  chemical  reaction  taking  place  and  the  transport  processes  occur¬ 
ring  in  the  system.  They  are  often  difficult  to  observe,  however,  as  the  experimental 
conditions  are  harsh  and  the  molecules  are  often  very  reactive, for  example,  free  radicals. 
Here,  we  present  an  ab  initio  technique,  which  is  within  the  Born-Oppenheimer  frame¬ 
work,  for  predicting  spectra  of  small  molecules  as  a  function  of  temperature.  It  uses  a 
combination  of  the  precise  and  effiecient  Complete  Active  Space  Self  Consistent  Field 
(CASSCF)  method  for  determining  electronic  wavefunctions,  with  a  fast  time-dependent 
ware  packet  technique  to  account  for  the  effects  of  nuclear  motion.  These  techniques 
are  generally  applicable  to  a  large  range  of  molecules,  and  may  be  applied  to  calculate 
spectra  of  samples  at  any  temperature.  Vibrational  and  electronic  spectra  of  C2H  are 
calculated  as  an  example  application  of  this  technique. 

The  ethynyl  radical  »C=  CH  is  believed  to  be  an  important  species.  Unassigned 
spectral  lines  in  the  interstellar  medium  1-4  have  been  attributed  to  C2H;  it  is  thought  to 
be  abundant  in  carbon  rich  stars;5'8  and  to  be  important  as  a  chemical  reaction  inter¬ 
mediate,9’  10  especially  on  Jupiter.11  In  the  laboratory  it  is  very  difficult  to  produce  and 
positively  identify  C2H.  Common  methods  of  production  include  photolysis12  of  C2H2 
and  electric  discharges13,  14  either  in  C2H2  or  over  polyacetylene.15  Recently,  a  summary 
of  production  methods  has  been  given  by  Laufer.7  The  most  important  signatures  of  the 
molecule  are  its  four  rotational  lines  near  87.3  GHz1, 12, 16, 17  and  its  CC  str*  tching  fun¬ 
damental13,  18  at  1848cm-1.  Jacox13  observed  the  CH  stretching  fundamental  for  C2H 
trapped  in  an  inert  matrix  at  14K  to  be  at  3612cm-1.  Under  similar  conditions  Graham, 
Dismuke  and  Weltner12  observed  the  lowest  lying  electronic  transition,  A(2n)<-  X(2£+), 
to  be  a  poorly  resolved  band  in  the  red  near  10000cm-1.  This  band  is  strongly  solvent 
dependent  and  in  the  gas  phase  is  both  observed15  and  predicted19  to  have  its  origin 
transition  near  3600cm-1.  Carrick,  Merer  and  Curl15  have  taken  high  resolution  spectra 
in  the  3600-4200cm-1  region  and  have  identified  five  vibronic  levels.  None  of  these  lev¬ 
els  were  identified  as  the  ground  state  CH  stretching  fundamental  i/,  which  was 
presumed  to  be  swamped  in  intensity  by  the  electronic  transition.  Interpretation  of  the 
high  resolution  spectra  is  very  difficult  as  there  is  strong  Renner-Teller  coupling  in  the  A 
state,  strong20  coupling  between  the  A  and  X  states  once  the  molecule  bends,  and  strong 
sp  v.bit  coupling.  ; 
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Shih.  Peyerimhoff  and  Buenker  have  made  ab  initio  calculations  of  both  the  verti¬ 
cal  electronic  spectrum19  and  the  emission  spectrum21  of  C2H.  They  predicted  that  the 
origin  and  center  of  the  A«-  X  absorption  band  should  be  around  4000cm-1  and 
5800cm'1,  respectively.  Also,  they  predicted  that  the  band  should  show  a  progression  in 
the  excited  state  CC  stretch  vibration  i/3.  The  majority  of  the  excited  states  were 
predicted  to  have  long  CC  bond  lengths  and  bent  equilibrium  geometries,  accounting  for 
the  vast  differences  observed  between  absorption  and  emission  spectra.10,  21 

The  vertical  electronic  spectral  calculations  of  Shih,  Peyerimhoff  and  Buenker19 
predict,  at  OK,  the  total  intensity  and  band  center  of  the  absorption.  They  give  only 
qualitative  information  on  the  bandwidth  and  no  information  on  either  structure  within 
the  band  or  on  the  change  of  band  shape  with  temperature.  Experimentally,  only  a  few 
lines  in  the  low  frequency  tail  of  the  A<-  X  absorption  band  have  been  observed  in  the 
gas  phase.15  This  paper  describes  a  calculation  of  the  principal  intensity  distribution  of 
the  A <-  X  absorption  band  at  3000K.  This  is  shown  to  be  dramatically  different  from 
the  OK  spectrum.  High  temperature  information  is  important  to  combustion7, 10  studies, 
as  well  as  to  the  problem  of  the  heating  of  a  space  probe  entering  the  Julian  atmo¬ 
sphere.22,  23 

We  use  the  Complete  Active  Space  Self  Consistent  Field7,10,24,25  (CASSCF) 
method  to  calculate  the  adiabatic  potential  surfaces,  fitting  this  data  to  a  functional 
form  of  the  potential.  Infrared  and  electronic  spectra  are  then  calculated  using  the  time 
dependent  quantum  formalism  of  Reimers,  Wilson  and  Heller.26  In  order  to  apply  this 
formalism  we  make  the  assumption  that  the  full  A(2IT)<-  X(2£+)  absorption  spectrum 
can  be  written  as  the  sum  of  spectra  calculated  for  each  component  of  the  bent  A  state, 
A(2A')*-X  and  A(2A")—  X.  The  results  show,  at  low  resolution,  how  the  gross  features 
of  the  absorption  contour  change  shape  with  temperature.  The  main  changes  in  the 
intensity  distribution  with  increasing  temperature  are  seen  to  be  a  blue  shift  due  to  the 
very  high  bend  vibrational  frequency  of  the  A(2A')  state  and  a  broadening  due  to  the 
considerable  anharmonicity  of  the  potential  surfaces. 

The  CASSCF  method24,  25  is  chosen  to  generate  the  potential  surfaces  for  the  fol¬ 
lowing  reasons.  First,  the  method  employs  a  fixed  configuration  list  at  all  geometries, 
thus  being  capable  of  high  precision.  This  is  an  essential  consideration  when  the  poten¬ 
tials  are  to  be  fitted  to  a  functional  form.  Second,  the  method  leads  to  proper  dissocia¬ 
tion  and  tends  to  treat  the  surfaces  with  uniform  accuracy  even  for  large  displacements 
from  equilibrium.  Finally,  the  CASSCF  method  is  sufficiently  economical  to  permit  a 
large  number  of  energy  and  dipole  moment  evaluations  on  both  the  X  and  A  surfaces. 
More  elaborate  methods  of  including  correlation,  such  as  multi-reference  singles  plus 
doubles  configuration  interaction,  would  require  a  very  large  number  of  references  to 
describe  the  entire  surfaces  equivalently.  Selection  methods,  which  would  make  the  cal¬ 
culation  feasible,  tend  to  reduce  significantly  the  precision  of  the  surface. 

Traditionally,  spectra  would  be  calculated  by  first  solving  the  Watson  Hamil¬ 
tonian,27  assuming  adiabatic  potential  surfaces  and  no  spin-orbit  coupling.  The  extra 
couplings  would  then  be  included  in  a  variational  calculation  and  the  observed  energies 
and  intensities  calculated.  This  is  an  extremely  ambitious  calculation,  as  many  hundreds 
of  vibronic  states  contribute  to  the  300QK  spectrum,  and  each  of  these  states  will  have 
complicated  PQR  rotational  structures.  At  the  moment  such  a  calculation  is  not  feasible. 
Consider,  by  way  of  example,  the  well-studied  case  of  the  water  molecule.  Compared  to 
C2H,  it  has  relatively  low  anharmonicity,  yet,  using  conventional  normal  mode 
approaches,  it  is  currently  only  possible  to  determine  the  lowest  20  vibrational  levels  of 
water  accurately28  (though  use  of  local  mode  basis  sets  significantly  improve  this  calcula¬ 
tion29  ). 

No  vibrational  eigenfunctions  are  required  in  the  approach  of  Reimers,  Wilson  and 
Heller.26  Here  the  spectrum  is  given  as  the  Fourier  transform  of  the  time  correlation 
function  obtained  by  evaluating  the  phase  space  average  of  the  overlap  of  two  gaussian 


wave  packets.  One  of  these  wave  packets  propagates  on  the  ground  potential  energy 
surface  while  the  other  propagates  on  the  excited  potential  energy  surface.  The  frozen 
gaussian  approximation26,  30  to  wave  packet  dynamics  is  used  in  all  calculations,  and  a 
classical  approximation26  to  the  temperature  operator  is  used  in  all  calculations  at 
nonzero  temperature.  Such  a  calculation  is  termed  FC:  here  F  stands  for  "Frozen"  and  C 
stands  for  ''Classical".  In  view  of  all  of  the  other  approximations  used  in  this  calculation, 
it  is  felt  that  improving  upon  either  of  the  FC  approximations  is  not  justified  at  this 
time. 

All  of  the  approximations  made  are  insignificant  for  band  contour  calculations  as 
the  contour  is  determined  by  only  the  first  few  femtoseconds  of  the  time  correlation 
function.  Band  contour  calculations  are  also  very  fast,  typically  requiring  less  than  a 
minute  of  DEC  VAX-780  computer  time.  The  spectra  reported  here  are  at  60cm-1  reso¬ 
lution  corresponding  to  a  propagation  time  of  530fs  Non-adiabatic  effects,  spin-orbit 
coupling,  wave  packet  deformations,  etc,  all  affect  the  dynamics  at  long  times,  causing 
the  calculated  lines  to  be  split  and  the  intensity  pattern  to  broaden.  Clearly,  60cm-1 
resolution  is  too  narrow  (the  spin-  orbit  coupling  alone  is  known  to  be  52cm-1,  for  exam¬ 
ple),  but  it  is  best  to  present  spectra  which  are  known  to  be  too  sharp  rather  than  spec¬ 
tra  which  are  thought  to  be  too  broad.  A  typical  calculation  at  this  resolution  takes 
0.5-2.0  hours  of  VAX-780  computer  time. 

II.  CALCULATION  METHODS 

A.  Ah  Initio  Calculations 

All  of  the  X,  A(2A)  and  A(2A")  surfaces  are  calculated  over  a  large  geometry  range 
using  the  CASSCF  method.  The  X  state  is  defined  over  the  range  accessible  thermally 
at  a  temperature  of  3000K,  while  the  A(2A')  and  A(2A")  surfaces  are  defined  over  the 
entire  energy  range  accessible  to  trajectories  starting  vertically  above  this  ground  state 
distribution. 

The  CASSCF  potential  energy  surfaces  are  generated  using  a  double-zeta  plus 
polarization  (DZP)  basis  set  of  cartesian  gaussian  functions.  For  carbon  we  use  the 
[4«  2 p  j  contraction31  of  the  (9s  5 p  j  Huzinaga  set,32  augmented  by  a  d  function  with  an 
orbital  exponent  of  0.75.  For  hydrogen,  we  use  the  [3s]  contraction33  of  the  |5s]  Huzi¬ 
naga32  set,  multiplied  by  a  scale  factor  of  1.2  and  augmented  by  a  p  function  with  an 
orbital  exponent  of  1.728. 

In  our  CASSCF  calculation  the  seven  outermost  valence  electrons  are  correlated  in 
seven  active  orbitals  -  five  of  a  symmetry  and  two  of  a"  symmetry.  All  calculations 
(even  at  C2v  geometries)  are  carried  out  in  C,  symmetry  to  ensure  that  no  discontinuities 
arise  on  the  A  surfaces  upon  bending.  Making  all  possible  arrangements  of  the  seven 
electrons  in  the  seven  orbitals  consistent  with  spin  and  symmetry  restrictions  results  in 
404  configurations  for  the  A(2A')  state  and  380  configurati  >ns  for  the  A(2A")  state.  Cal¬ 
culations  with  both  smaller  and  larger  active  spaces  indicates  that  these  CASSCF  calcu¬ 
lations  are  reasonably  well  converged. 

The  CASSCF  dipole  moment  surface  is  used  to  compute  the  infrared  spectrum  of 
the  X(2E+)  ground  state.  Analogous  CASSCF  calculations  on  the  ozone  molecule34  pro¬ 
duce  intensities  for  the  fundamental  and  lowest  combination  and  overtone  bands  that 
are  in  excellent  agreement  with  experiment.  Further  improvement  in  the  dipole  moment 
function  would  require  the  inclusion  of  a  diffuse  d  function  in  the  basis  set,  as  well  as  a 
more  extensive  treatment  of  the  electron  correlation. 

Since  the  CASSCF  wavefunctions  for  the  X  and  A  states  are  defined  in  different 
orbital  bases,  we  can  not  use  these  wavefunctions  to  deduce  the  electronic  transition 
moment  functions.  Instead  we  use  transition  moments  based  upon  Cl  wavefunctions 
defined  in  terms  of  the  ground  state  Hartree-Fock  orbitals.  The  transition  moments  are 
rather  insensitive  to  electron  correlation  owing  to  the  fact  that  the  dominant 


configurations  of  the  X  and  A  states  differ  by  one  spin  orbital.  Hence,  the  transition 
moments  are  calculated  using  only  small  Cl  wavefunctions  with  selection  based  on  the 
magnitude  of  the  diagonal  hamiltonian  matrix  elements.  The  calculated  oscillator 
strength  at  the  equilibrium  geometry  is  in  excellent  agreement  with  more  extensive  sin¬ 
gles  plus  doubles  Cl  treatments  as  well  as  with  the  results  of  Shih  Peyerimhoff  and 
Buenker.19  Note  that  in  this  calculation  only  three  cuts  through  the  full  transition 
moment  profile  are  determined. 

B.  Spectroscopic  Analysis 

The  spectroscopic  programs  are  very  general  and  can  be  written  in  a  variety  of  co¬ 
ordinate  systems.  Typically  they  are  written  in  cartesian  co-ordinates.  These  are  concep¬ 
tually  the  simplest  and  require  no  separation  of  rotational  and  vibrational  variables.  For 
proteins,  liquids  and  moderate  sized  molecules  they  are  clearly  the  best  possible  co¬ 
ordinates  but,  for  triatomic  molecules,  they  have  the  disadvantage  of  containing  three 
times  the  number  of  variables  that  a  purely  vibrational  set  of  co-ordinates  would  con¬ 
tain.  Initially  we  planned  to  run  thawed  gaussians,26-  35  for  which  this  problem  is 
accentuated  as  the  second  derivatives  of  the  potential  surface  are  then  required.  Also 
thawed  gaussians  require  special  treatment  to  include  rotational  degrees  of  freedom.36,  37 
Finally,  cartesian  co-ordinate  potential  functions  require  a  large  number  of  square  root 
and  trigonometric  function  evaluations,  which  significantly  slow  the  programs. 

A  common  alternative  to  cartesian  co-ordinates  is  valence  co-ordinates  (bond 
lengths  and  angles).  Their  use  introduces  the  rigid  rotor  approximation.  The  potential 
energy  then  takes  a  much  simpler  form,  at  the  expense  of  complicating  the  kinetic  energy 
operator.  Another  approach  is  to  use  normal  co-ordinates.  They  again  restore  the 
kinetic  energy  operator  to  a  simple  form  but  are  usually  regarded  as  poor  variables  in 
which  to  expand  a  potential  surface.  Here  a  simple  non-linear  transformation  is  used  to 
construct,  from  the  normal  co-ordinates,  a  new  set  of  co-ordinates  which  closely  resemble 
the  valence  co-ordinates.  The  spectroscopic  programs  can  thus  be  written  entirely  in 
normal  co-ordinates.  In  this  approach  a  large  part  of  the  anharmonicity  of  the  spectra  is 
attributed  to  the  intrinsic  curvilinear  nature  of  the  molecular  potential  surface.  This  is 
equivalent  to  the  unphysical  concept  of  "kinetic  energy  anharmonicity",  which  is  com¬ 
monly  referred  to  in  the  literature  of  calculations  that  use  valence  co-ordinates. 

Using  normal  co-ordinates,  the  basic  formalism26  for  the  infrared  or  electronic  cross 
section  a  T  (w),  at  temperature  T ,  becomes  the  Fourier  transform 

=  “jp-  (l-e-pTu)  J_l  CT  (t  )eiwl  dt  (2.1) 


of  the  time  correlation  function 
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Cr(0  =  Q  1  f_godnq0dHp0  (xqj.J  e  e  * 1  ‘  % '  e  ,Nt,lrp  \  xqoPo)  (2.2) 

where  /3  =  (JfcT  )-1,  k  is  Boltzmann’s  constant,  c  is  the  speed  of  light,  Q  is  the  quantum 
partition  function, 


Q  =  Xoo<inq0<iBp0(X0»l>ol  e'0*'  1 


(2.3) 


|  APoq0)  *s  a  gaussian  wave  packet  initially  centered  at  position  q0  and  momentum  p0,  n 
is  the  dipole  or  transition  moment  operator,  and  //,  and  Hf  are  the  initial  (X)  and  final 
(A)  state  hamiltonians,  respectively.  Here  we  take  q0  as  a  vector  of  dimensionless  nor¬ 
mal  co-ordinates  and  p0  as  their  conjugate  momenta.  The  classical  temperature  approxi¬ 
mation26  replaces  the  operator  e  '  by  the  number  e  fiE'  where  E ,  is  the  classical 


energy  of  the  trajectory.  This  approximation  arises  by  making  a  Taylor  expansion,  in 
powers  of  /3,  of  the  effect  of  the  quantum  operator  on  the  wavefunction  and  then  retain¬ 
ing  only  the  leading  term.  It  reduces  the  full  partition  function  to  the  classical  partition 
function, 


Q  =  f_00d”ti od"Po  e  fiE‘ 

(2.4) 

and  allows  the  integrand  in  equation  (2.2)  to  be  written 
e  ($,■  |  $  f  ^  where 

as  Boltzmann  weighted  overlap 

|  4*,)=  |  xqaPo) 

(2.5) 

and 

1  */)  =  e~'W/</r^  1  Xq<)p0)  ‘ 

(2.6) 

Equation  (2.5)  is  interpreted  as  propagating  the  wave  packet  on  the  initial  (X)  potential 
surface  to  time  t  and  then  then  operating  upon  the  wave  function  with  the  dipole 
moment  operator,  while  equation  (2.6)  is  interpreted  as  propagating,  on  the  final  (A) 
state,  the  result  of  the  dipole  moment  operating  on  the  wave  packet. 

In  this  paper  we  use  the  frozen  gaussian  approximation30  to  propagate  the  wave 
packets.  This  constrains  the  width  of  the  gaussian  function  to  be  constant  in  time.  The 
center  of  the  gaussian,  initially  p0  and  q0,  then  moves  according  to  Hamilton's  equations 
of  motion,  and  the  wave  packet’s  phase  contains  a  term  reflecting  the  classical  action  as 
well  as  a  zero  point  motion  correction  term.  It  is  thus  no  more  difficult  to  implement 
than  is  classical  molecular  dynamics  and  yet  is  an  approximation  which  retains  most  of 
the  essential  features  of  quantum  mechanics.  Tunneling  can  not  be  directly  accounted 
for  as  this  demands  that  the  wave  packet  splits  into  pieces.  Features  that  are  retained 
by  this  approach  include  barrier  penetration,  the  Heisenberg  uncertainty  principle,  zero 
point  motion,  interference  effects  and  basic  operator  commutation  relationships. 

For  linear  triatomic  molecules  such  as  C2H,  there  are  eight  phase  space  (position 
and  •uomentum)  variables  to  integrate  over  in  equations  (2.2)-(2.3).  These  integrals  are 
evaluated  numerically  using  a  product  Gauss-Hermite  quadrature.38  This  quadrature  is 
accelerated  by  only  running  trajectories  with  significant  weight  and  by  transforming  to 
variables  which  better  reflect  the  X  potential  surface. 

The  high  temperature  formalism  presented  here  requires  no  adaptation  to  calculate 
spectra  at  OK.  At  this  temperature  the  classical  phase  space  available  is  just  one  point, 
and  the  gaussian  wave  packet  centered  on  this  point  is  just  the  harmonic  approximation 
to  the  ground  state  wavefunction.  This  formalism  thus  reduces  to  the  eigenstate  formal¬ 
ism  of  low  temperature  wave  packet  d /namics.39  The  classical  approximation  is  thus 
accurate  at  both  low  temperatures  ( kT  /JTw  <  0.5)  and  high  temperatures  (  kT  /iTu  > 
5).  At  3000K,  the  vibrations  of  C2H  fall  into  the  intermediate  region  with  kT  ftTu  = 
0.6,  1.1  and  4.2  for  the  bend  (t/2),  CC  stretch  (v 3)  and  CH  stretch  (i/j)  vibrations, 
respectively. 

m.  FITTING  THE  AB  INITIO  DATA 

A.  Functional  Form  of  the  Potentials 

The  CASSCF  data  for  each  state  S  =  X,  A(2A)  or  A(2A")  are  fitted  to  two  poten¬ 
tial  functions.  Both  of  these  functions  have  the  same  form  but  they  are  written  as  func¬ 
tions  of  different  position  variables.  One  function,  Fjsl(Ar),  is  a  function  of 
AT  =  {Arc*  ,A0,A rcc  },  the  change  in  the  valence  co-ordinates  T  =  {rclj  ,0,rcc  }  from 
the  X  state  equilibrium  geometry  r*x'  ,rjx^},  where  rCH  and  rcc  are  the  CH 


and  CC  interatomic  distances,  respectively,  and  6  and  if  are  the  CCH  bond  angle  and 
its  equilibrium  value,  respectively.  The  other  function.  V^JS|(q).  is  a  function  of  q,  a  set 
of  dimensionless  norma]  co-ordinates  of  the  X  state,  q  =  {9],?2.?j}  where  qx  is  predom¬ 
inately  the  CH  stretch,  q2  is  the  bend  and  q2  is  predominately  the  CC  stretch.  Since 
C2H  is  linear,  the  bend  co-ordinates  A0  and  q2  have  two  degenerate  components.  If  <t>  is 
an  angle  which  specifies  the  angle  between  the  plane  of  a  bent  molecule  and  the  plane 
perpendicular  to  the  original  molecular  axis  then  we  take  these  components  as  the  carte¬ 
sian  co-ordinates 

x  =  A0,  =  sinA0  cos<p  and  y  =  A0t  -  sinA0sin 4>  (3.1) 

and  the  normal  co-ordinates 

=  9 2cos<t>  and  q2k  =  q2  sin<£  .  (3.2) 

Note  that  the  potentials  are  cylindrically  symmetric  and  independent  of  <j>.  The  poten¬ 
tials  Vjsl(Ar)  are  written  as  functions  of  A0  rather  than  sinA0  as  A 0  has  superior  boun¬ 
dary  properties  at  large  6.  Only  Vjsi(q)  are  used  in  trajectory  calculations.  The 
K„^(Ar)  surfaces  are  presented  for  comparison. 

The  valence  co-ordinate  potentials  Vjsl(Ar)  are  written  in  a  form  containing  27 
free  parameters.  Removing  the  superscript  [S]  for  clarity,  this  form  is  written  as 

Vj  ^(Ar)  ~  V0  +  DHHxk  +  Dhc  xk  xe  +  Dccxc  +  DB{xk-\)  +  EB  (xk  —\^  ^-3^ 


where 

xk  =  1  -  exp(-  aH  6h  )  ,  (3.4) 

xc  =  1  -  exp(-ac<5c)  ,  (3.5) 

xh  =  exp(a„  b2)  ,  (3.6) 

a„  =  a0"  +  aUh  +  a"c  +  aBBb2  ,  (3.7) 

°c  =  a0c  *r  qh  h  +  accc  -f  aBB  b 2  ,  (3.8) 

aB  =  a0“  ■+  a*h  +  a‘c  -f  aBBb 2  ,  (3.9) 

&h  ~  h  +  Dhh  h  ■+■  Dhc  c  +  8  BB  6 2  ,  (3.10) 

8C  =  c  +6„ch  +  Dccc  ■+  8bb  b 2  ,  (3.11) 

E„  =  Ef  +  E»c  +  E»c2  ,  (3.12) 

and  the  displacements  from  the  equilibrium  geometry  of  state  5  are  given  by 

h  =  ArCH  -  (r^-rj*1)  ,  (3.13) 

c  =  A rcc  -  (rjf-rj?1) ,  (3.14) 

and 

6 2  =  (A0)2  (3.15) 


where  Tc$  and  tJP  are  the  equilibrium  bond  lengths  of  state  S . 

The  stretch  potentials  are  basically  Morse  functions.  Due  to  symmetry,  the  bend 
potential  must  be  an  even  function  of  A#  and  we  find  it  convenient  to  write  this  as  a 
function  of  (A/0)2.  The  basic  form  used  for  the  bend  potential  is  a  combination  of  a 
gaussian  and  a  rising  exponential.  It  is  flexible  enough  to  allow  the  equilibrium 
geometry  to  be  non-linear,  if  necessary,  and  fits  the  C2H  surfaces  much  better  than  does 
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potential  surfaced  based  upon  Taylor  expansions.  Other  terms  are  added  to  these  basic 
functions  to  allow  them  to  change  curvature  slightly  at  locations  away  from  equilibrium. 


Many  other  functional  forms  of  this  potential  have  been  tested  and  found  to  be 
unsatisfactory.  Taylor  expansion  type  potentials,  including  the  potentials  of  Carney, 
Curtiss  and  Langhoff28  (an  expansion  in  the  Simon-Parr-Finlan  coordinate40  Ar  /r)  and 
a  generalized  Hoffaker  potential41  (an  expansion  in  the  Morse  variable  ]-e_oAr )  are  too 
flexible.  They  are  difficult  to  control  as  their  basic  shape  can  change  rapidly  over  short 
distances  producing  nasty  saddle  points  at  relatively  low  energies.  The  function  used 
here  is  not  as  flexible  as  these  forms  and  its  saddle  points  tend  to  be  at  high  energies  and 
outside  the  region  of  the  potential  being  fitted.  It  still  has  sufficient  flexibility  to  inter¬ 
polate  the  potential  surface  accurately. 

Valence  co-ordinates  reflect  the  shape  of  the  potential  surface  much  better  than  the 
normal  co-ordinates,  which  are  just  the  tangential  components  of  the  valence  co¬ 
ordinates  at  the  equilibrium  geometry.  The  transformation  from  normal  to  valence  co¬ 
ordinates  has  been  specified  up  to  third  order  by  Hoy,  Mills  and  Strey.42  Our  approach 
is  to  include  only  the  linear  and  dominant  non-linear  terms  thus  generating  from  the 
norma)  co-ordinates  a  new  set  of  co-ordinates  s  =  {sH  ,sB  ,sc  }  (where  again  sBA  =  sB  cos</> 
and  sBB  =  sB  sin^  )  that  closely  resemble  the  valence  co-ordinates  Ar.  These  are  also  a 
good  set  of  variables  in  which  to  expand  the  potential,  and  are  defined  as: 


sh  ~  Lhh  0i  +  L  HC  9  3  +  ^HBB  9  2  “t"  L HBBBB  9  2  » 
SC  ~  Lch  0i  +  L  cc  03  +  Lcbb  q 2  +  LCBBBb  0 2  > 
sb  2  ~  Lbb  20j2  +  Lbbbbb  q 2  • 


(3.16) 

(3.17) 

(3.18) 


The  normal  co-ordinate  potentials  VnS|(q)  are  defined  in  an  identical  fashion  to  the 
valence  co-ordinate  potential,  except  that  s„  ,  sB2  and  sc  replace  A TCH  ,  (A0)2  and  A Tcc , 
respectively,  in  equations  (3.13)-(3.15). 

In  the  equations  (3.16)-(3.18)  the  five  non-linear  L  terms  are  treated  as  free  param¬ 
eters  and  are  refined  along  with  the  potential  constants.  The  five  linear  L  terms  define 
the  direction  and  length  of  the  dimensionless  normal  co-ordinates  q.  If  q  are  to  be  the 
normal  co-ordinates  of  Vjxl(q),  then  the  linear  terms  are  dependent  upon  the  other  force 
constants. 


B 

in 


.  Fitting  the  Potential  Co-efficients 

The  potential  coefficients  are  obtained  by  refining  their  values  to  minimize  the  error 
the  fitted  potential  Vnr  according  to 

,,2  _  (  Vf/r  -  Vcasscf  )2 


Vcass  ~r  4  10cm  1 


(3.19) 


where  the  sum  is  over  the  ab  initio  CASSCF  data  points.  This  form  of  the  error  func¬ 
tion  is  superior  to  a  least  squares  form  as  it  does  not  over  emphasize  the  high  energy 
points.  For  the  X  state  potential  energy  surfaces  all  of  the  force  constants  are  adjusted 
simultaneously  to  fit  all  of  the  CASSCF  data.  For  the  A  state  surfaces  a  different 
approach  is  taken  and  the  stretch-only  force  constants  are  fitted  to  the  linear  data  points 
only.  All  of  the  A  surfaces  thus  have  the  same  stretching  potential.  There  are  four  dif¬ 
ferent  bending  potentials:  for  each  of  the  two  components  of  the  A  state  there  are  two 
functions  generated,  one  in  the  valence  co-ordinates  and  the  other  in  the  normal  co¬ 
ordinates. 

In  Table  I  the  potential  constants  are  given  and  in  Table  II  the  accuracy  of  the  fits 
are  summarized.  Table  II  gives,  for  each  of  the  seven  fits  made,  the  total  number  of 
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E  I.  Potential  Coefficients.  ^ 


Constant 

VJxl(q) 

V’JA(,A')i(q) 

VJA(lA")l(q) 

V,w(Ar) 

vJa(2a")!(At) 

*ch  /A 

1.056180 

1.059800 

1.059800 

1.056180 

1.059800 

1.059800 

Tcc  /A 

1.232760 

1.314430 

1.314430 

1.232760 

1.314430 

1.314430 

/A  ^ 

-0.148763 

-0.380053 

0.068499 

'HBBBB  /A 

-0.049642 

-0.052920 

-0.330783 

CBB  /A 

-0.004142 

-0.005807 

-0.008128 

‘CBBBB  /A 

0.001944 

-0.002816 

-0.003152 

Lbbbbb 

0.011084 

-0.003409 

0.005216 

ao  /A  -  ’ 

1.995370 

1.790790 

1.790790 

1.737230 

1.790790 

1.790790 

o.*/A  "2 

0.383248 

-0.001382 

-0.001382 

-0.047890 

-0.001382 

-0.001382 

®c7A  2 

-0.292976 

0.190334 

0.190334 

-0.045747 

0.190334 

0.190334 

a/s/A  1 

-0.906130 

-0.746859 

-0.755775 

-0.071663 

-0.015817 

-0.053669 

«o7A 

1.848800 

1 .973590 

1.973590 

1.918900 

1.973590 

1.973590 

o„c/A  2 

0.023466 

-0.475799 

-0.475799 

0.016937 

-0.475799 

-0.475799 

Oce/A  2 

-0.191945 

-0.018958 

-0.018958 

-0.102989 

-0.018958 

-0.018958 

a  *8  /A  1 

0.841101 

-0.308327 

-0.155382 

-0.728090 

-0.441608 

0.136459 

ao* 

0.289789 

1.867800 

0.165656 

0.274980 

2.074430 

0.702795 

°2 

-0.009747 

-1.091110 

-0.316241 

0.042571 

-1.886850 

-0.042457 

°c 

1.323840 

2.139090 

-1.191640 

-0.378518 

4.616360 

-1.145650 

< 

0.046013 

0.469322 

0.557571 

0.022000 

-0.044138 

-0.110016 

Ohm 

-0.040104 

-0.002121 

-0.002121 

0.009888 

-0.002121 

-0.002121 

Due 

0.137949 

-0.020134 

-0.020134 

-0.026061 

-0.020134 

-0.020134 

Kb!  A 

-0.179276 

0.119968 

-0.381256 

0.018716 

-0.008092 

0.010861 

-0.001134 

0.050000 

0.050000 

-0.027138 

0.050000 

0.050000 

Dec 

0.024080 

0.003973 

0.003973 

0.003405 

0.003973 

0.003973 

f>BB  /  A 

-0.053271 

0.104576 

0.201632 

0.068520 

0.116599 

-0.033598 

Vr0  cm  ' 

0 

3600 

3600 

0 

3600 

3600 

P««  /cm  1 

41671.5 

54287.8 

54287.8 

60039.8 

54287.8 

54287.8 

D„c  /cm  1 

11344.8 

3228.6 

3228.6 

-11482.3 

3228.6 

3228.6 

c  /cm  *A  1 

114035.0 

66573.4 

66573.4 

100870.0 

66573.4 

66573.4 

,  /cm_,A  -1 

11056.6 

153.6 

8272.9 

10591.4 

75.0 

5341.3 

E f  /cm  1 

248.1 

-7386.5 

-9076.4 

0.0 

-6961.1 

0.0 

Ef  /cm1 

362784.0 

-62222.6 

49633.0 

5227.4 

-46671.6 

43236.5 

Ecc  /cm  1 

-486440.0 

53941.5 

-20208.9 

-409647.0 

29038.2 

7901.8 

lata  points,  the  number  of  adjustable  parameters  n ,  the  final  value  of  \2,  the  average 
ibsolute  error  as  a  percentage  of  the  average  energy,  and  the  largest  recorded  error, 
tlore  detail  and  the  CASSCF  data  points  used  are  available  on  request  from  the 
luthors.  Most  of  the  data  points  are  within  10000cm-1  of  the  minima  of  the  potential 
surfaces,  though  there  are  some  points  on  each  surface  in  excess  of  15000cm"1  energy. 
CJeometries  with  CCH  angles  as  low  as  100°  are  included  in  this  data.  Both  forms  of 
:he  potential  surface  fit  the  data  well.  The  distribution  of  the  errors  in  the  fits  are  ran- 
Jom  and  there  is  little  correlation  between  corresponding  normal  and  valence  co-ordinate 
surface  errors.  Typically  the  functions  fit  the  data  to  a  tolerance  of  0.5%,  with  the 
greatest  error  observed  for  any  one  point  being  331cm-1. 

All  energies  reported  in  Table  1  and  elsewhere  are  with  respect  to  the  minimum 
energy  of  the  X  state  which  is  calculated  to  be  -76.2636  Hartrees.  The  calculated  spac¬ 
ing  between  the  minima  of  the  X  and  A  states,  ,  is  5835cm-1,  in  contrast  to  pure 
SCF  calculations  which  actually  place  the  A  state  at  a  little  lower  energy  than  the  X 
state.  Experimental  results15  place  VjAl  near  3600cm-1  in  reasonable  agreement  with 
multi-reference  singles  plus  double  Cl  results.19  Apparently  the  CASSCF  treatment  with 
five  active  a'  orbitals  and  two  active  a"  orbital  provides  a  somewhat  better  description  of 
the  X  state  than  of  the  A  state.  We  have  computed  V(jA'  using  Cl  wavefunctions 
(including  all  singles  and  doubles  replacements  from  the  Hartree-Fock  reference)  in  both 
our  DZP  basis  and  an  extended  van  Duijneveldt  gaussian  basis43  (|13s7p2d]  contracted 
to  [7s4p2rfj).  At  the  Cl  singles  plus  double  level,  PjAl  is  2173cm"1  and  2462cm  1  for 
the  DZP  and  extended  bases,  respectively.  These  energies  increase  substantially  to 

2882cm"1  and  3046cm-1  when  a  correction  is  applied  for  quadrupole  excitations.44 

Hence,  as  more  of  the  valence  correlation  energy  is  recovered,  a  larger  proportion  of  the 
differential  correlation  is  recovered  as  well.  It  is  likely  that  a  fully  converged  Cl  calcula¬ 
tion  will  produce  a  VjAl  in  good  agreement  with  the  experimental  observation  of  near 
3600cm"1,  so  we  use  this  value  in  our  calculations  instead  of  5835cm-1. 

From  the  valence  co-ordinate  potential  V'Jx*(Al*)  the  dimensionless  normal  co¬ 
ordinates  q  are  deduced.  These  co-ordinates  have  the  frequencies 

i/j  =  3623cm  \  =  486cm  1  and  i/3  -  1971cm  where  =  u;,  /2tr  c  .  Dimensionless 

normal  co-ordinates  are  produced  in  the  usual  fashion  by  weighting  the  dimensioned  nor¬ 
mal  co-ordinates  by  (TTw,)  .  The  direction  and  length  of  these  co-ordinates  are  then 
obtained45  from  the  eigenvectors  of  the  matrix  of  the  mass  weighted  second  derivatives 
of  Pjxl(Ar),  expressed  in  cartesian  co-ordinates.  This  information  is  contained  within 
the  parameters  LHH  =  0.099548A  ,  L„c  =  0.013497 A  ,  LCH  =  -0.011555A  , 
Lcc  =  0.051046A  and  LBB  =  0.288801.  All  of  the  normal  co-ordinate  potentials  are 
written  in  terms  of  these  dimensionless  normal  co-ordinates.  Note  that  as  all  of  the 
potential  constants  are  refitted  to  the  data  when  Vjs*(q)  are  generated,  the  normal  co¬ 
ordinate  potentials  need  not  have  the  same  vibration  frequencies  as  Pe^(Ar). 

C.  Dipole  and  Transition  Moments 

Equation  (2.6)  requires  the  result  of  the  dipole  moment  operating  on  the  gaussian 
wave  packet  to  be  propagated  on  the  excited  potential  surface.  This  can  be  propagated 
without  further  approximation26,  46  if  /x  can  be  expanded  as  a  power  series  in  the  dimen¬ 
sionless  normal  co-ordinates  q.  In  principle,  all  terms  in  the  expansion  can  be  accommo¬ 
dated  but  at  present  only  constant  and  linear  terms  are  included  in  the  spectral  pro¬ 
grams.  When  higher  order  terms  are  small  in  comparison  to  these  terms  (this  does  not 
apply  to  the  CC  stretch  component  of  the  X  state  dipole  moment),  their  effect  may  be 
approximated  by  locally  expanding  p  as  a  linear  function. 

Vibrational  transitions  to  levels  with  E+  symmetry  are  polarized  parallel  (||)  to  the 
molecular  axis.  The  associated  dipole  moment,  //, ,  is  an  even  function  of  the  bend  coor¬ 
dinates.  Conversely,  transitions  to  levels  with  IT  symmetry  are  polarized  perpendicular 
(-*.)  to  the  molecular  axis  and  the  dipole  moment  fxt  is  an  odd  function  of  the  bend  co- 


TABLE  III.  Transition  Moment  and  Dipole  Moment  Coefficients,  in  Debyes 


|  Constant 

X 

A(2A')<-  X 

A(2A")<-  X 

A4  o 

0.725579 

0.549164 

0.549164 

1  Vh 

0.089610 

0.008069 

0.008069 

Vc 

0.026924 

-0.024602 

-0.024602 

!  A*  HH 

-0.001546 

-0.000791 

-0.000791 

j  A4  HC 

0.007610 

0.000451 

0.000451 

i  A*  cc 

0.000474 

-0.001439 

-0.001439 

A4  as 

0.014547 

-0.032514 

0.007478 

•  A'  HBB 

-0.002859 

0 

0 

|  A1  CBS 

0.023797 

0 

0 

i  A4  HHBB 

0.000146 

0 

A4  HCBB 

0.001957 

0 

o 

|  A*  CCBB 

-0.001607 

0 

0 

|  A'  BBBB 

-0.002209 

0.001615 

0.000713 

1 

0.168037 

0.319952 

1  A'hb 

0.006908 

0 

| 

1  A^cs 

-0.030209 

0 

A4  HHB 

0.003773 

0 

A1  HCB 

-0.003326 

0 

j  A4  cca 

0.000769 

0 

A*  BBS 

0.005092 

0.000382 

A4  HHBB 

0.000198 

0 

A4  HHCB 

-0.000189 

0 

| 

A1 HCCB 

-0.000041 

0 

1 

1 

A4  CCCB 

0.000560 

0 

j 

;  A4  HBBB 

0.001286 

0 

j 

A4  CBBB 

0.000629 

0 
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C2H  X  VIBRATIONAL  ABSORPTION 


Figure  5.  Vibrational  Absorption  Spectrum  calculated  using  either  the  FC 
AC  methods  and  either  the  full  dipole  moment  or  only  its  linear  components. 


or 
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large  number  of  vibrational  states  which  contribute  to  the  spectrum,  in  contrast  to  the 
previous  room  temperature  calculation  of  the  water  molecule  infrared  and  Raman  spec¬ 
tra.51  We  thus  calculated  the  vibrational  spectrum  at  3000K  using  classical  linear 
response  theory,52  as  formulated  by  Berens  and  Wilson.53  In  this  approach  the  classical 
autocorrelation  (AC)  of  the  dipole  moment  is  computed  and  Fourier  transformed  to  give 
the  infrared  spectrum.  At  the  harmonic  level  this  approach  gives  the  correct  quantum 
intensities  for  the  fundamental  vibrations  at  non-zero  temperature.53  Non-linear  com¬ 
ponents  of  the  dipole  moment  are  easily  included  in  this  calculation.  It  has  the  disad¬ 
vantages,  however,  of  incorrect  detailed  balance  for  overtone  and  combination  bands,  of 
an  incorrect  represention  of  the  anharmonic  coupling,  and  it  can  not  be  applied  as  T  -*  0. 

Figure  5  gives  the  spectrum  calculated  usi.  ?  both  the  wave  packet  (FC)  and  classi¬ 
cal  autocorrelation  (AC)  methods  at  3000K.  The  results  for  the  linear  dipole  are 
remarkably  similar  for  seven  of  the  eight  observed  bands,  but  the  intensity  profile  of  i/j 
is  reduced  by  a  factor  of  four  in  the  AC  calculation.  These  calculations  have  been 
repeated  using  a  single  harmonic  oscillator  and  the  results  compared  to  exact  quantum 
calculations.26  The  full  quantum  and  AC  results  are  almost  identical  and  have  an  inten¬ 
sity  of  only  75%  of  the  FC  result.  This  discrepancy  is  due  to  the  reduced  temperature 
kT  /Tu,  which  for  this  vibration  is  0.6.  in  the  intermediate  region  where  the  classical 
temperature  approximation  is  least  appropriate.  This  intensity  is  also  three  times  larger 
than  the  intensity  calculated  using  the  AC  technique  on  the  full  potential.  In  Figure  5, 
most  of  the  CH  stretch  intensity  is  dispersed  by  the  anharmonicity,  appearing  in  both 
the  combination  and  overtone  bands  and  in  the  elevated  base  line.  Such  a  dispersion  is 
unrealistic  and  the  FC  calculation  gives  a  better  description  of  the  overtone  and  combi¬ 
nation  band  intensities. 

Classical  autocorrelation  function  calculations  may  also  be  performed  using  the  full 
dipole  moment  profile,  and  the  results  are  also  given  in  Figure  5.  The  most  striking 
feature  of  this  spectrum  is  the  five  fold  increase  in  intensity  of  the  CC  stretch,  uz. 
Crude  time-independent  calculations  using  harmonic  potentials  and  only  the  fic  and 
Pcbb  terms  from  the  dipole  moment  expansion  give  results  in  excellent  agreement  with 
this  result.  The  extra  intensity  arises  from  many  hot  bands  which  are  not  resolvable 
from  each  other  at  60cm'1  resolution.  High-resolution  high-temperature  spectra  are 
expected  to  be  very  complicated. 

B.  Electronic  Spectra 

Calculated  spectra  over  the  temperature  range  0-3000K  are  given  in  Figure  6  for 
the  A(2A")<-  X  transition.  The  results  are  easily  understood  in  terms  of  a  standard 
eigenstate  picture.  At  OK  we  see  a  typical  Frank-Condon  progression  based  on  the  an 
allowed  (perpendicular)  origin  at  3400cm'1  This  frequency  is  less  than  the  energy 
difference  between  the  potential  minima  (3600cm'1)  because  of  the  zero  point  energy 
decrease  in  the  A(2A")  state.  Normally  a  frozen  gaussian  calculation  would  not  correctly 
detect  this  change  (as  a  thawed  gaussian  calculation  would)  and  here  it  has  been  added 
in  an  ad  hoc  fashion  before  Fourier  transforming.  In  reduced  units  the  A(2A")  potential 
minimum  is  displaced  by  =  -0.15  and  q  3  =  1.58  from  the  X  potential  minimum.  The 
progression  in  uz  thus  has  some  length  and  the  line  spacing  seen  in  the  spectra  is 
1620cm'1,  very  close  to  the  harmonic  A  state  CC  stretch  frequency  of  1655cm'1.  Just  to 
the  blue  of  the  third  strong  peak  is  the  first  line  of  a  progression  corresponding  to  an 
excitation  of  one  quanta  in  the  CH  stretch  Most  of  this  intensity  comes  from  the  1:2 
Fermi  resonance  between  the  stretching  motions.  The  tiny  peaks  at  1900cm'1  and 
3800cm'1  are  artifacts  due  to  the  approximate  propagation  methods  used. 

As  the  temperature  is  raised  to  1000K,  hot  bands  in  the  CC  stretch  vibration 
appear.  A  transition  from  vibrational  level  i  to  vibrational  level  f  of  the  n  th  vibration 
is  labeled  as  /n  .  The  origin  of  this  hot  band  is  thus  labeled  as  pi,  and  the  progression 
based  upon  it  is  350cm'1  to  the  red  of  the  ’told"  lines  based  upon  the  origin  O  (q1^3). 
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Fifurt  4.  Vibrational  Absorption  Spectrum  calculated  using  the  F  ’  method 
and  either  all  or  part  of  the  linear  dipole  moment  derivative  terms. 
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It  is  also  possible  to  approximate  the  center  of  the  absorption  band  by  examining 
what  difference  contours  lie  above  the  most  probable  regions.  At  T  -  0  this  gap  is 
5630cm-1.  At  3000K  the  most  probable  region  for  the  A(2A" ')*-  X  transition  lies  on  simi¬ 
lar  energy  contours  and  so  the  absorption  band  is  expected  to  broaden  but  not  shift  in 
band  center.  Due  to  the  repulsion  of  the  A(2A')  surface,  the  energy  gap  in  the  probable 
140° -160°  range  of  bond  angle  is  3-4000cm-1  larger  than  it  is  at  the  equilibrium  confi¬ 
guration.  The  A(2A')»-  X  absorption  band  is  thus  expected  to  be  significantly  blue 
shifted  at  3000K.  This  line  of  argument  can  be  formalized  and  a  simple,  almost  classical, 
algorithm  developed  for  the  determination  of  the  absorption  band  contour.48 

V.  CALCULATED  SPECTRA 

A.  Vibrational  Spectra  of  the  X  State 

Figure  4  shows  the  infrared  vibrational  spectra  of  the  X  state  of  C2H  calculated 
using  the  Frozen  gaussian  -  Classical  temperature  (FC)  approximations  to  wave  packet 
dynamics.  Only  the  linear  components  of  dipole  moment  function  are  used  in  these  cal¬ 
culations.  The  spectrum  at  OK  is  precisely  the  spectrum  which  would  be  obtained  using 
harmonic  approximations  and  time  independent  quantum  mechanics.  Here  the  intensi¬ 
ties  are  proportional  to  the  square  of  the  dipole  moment  derivatives.  There  is  a  large 
difference  between  the  calculated  and  observed  intensities.  From  the  low  temperature 
spectra  of  Jacox,13  the  CH  stretch  tq  is  seen  to  be  about  three  times  as  intense  as  the 
CC  stretch  vz,  while  Figure  4  shows  i/j  to  be  about  20  times  more  intense.  Also  the 
bend  absorption  i/2  is  calculated  to  be  quite  strong,  whereas  this  vibration  has  never 
been  directly  observed.  The  approach  taken  here  adequately  describes  the  low  tempera¬ 
ture  infrared  intensities  of  many  molecules,  including  ozone,34  but  seems  inappropriate 
for  C2H. 

The  discrepancy  between  the  intensity  ratios  for  vx  and  vz  is  due  largely  to  the 
neglect  of  non-linear  and  anharmonic  terms.  In  principle,  the  electrical  non-linearity 
may  be  included  directly  into  the  FC  calculations,26  but  it  is  much  more  difficult  to 
include  the  effects  of  mechanical  anharmonicity.46’ 49  Here,  we  choose  to  perform  an 
expensive  time-independent  calculation,  diagonalizing  the  vibrational  hamiltonian  matrix 
to  account  for  the  mechanical  anharmonicity,  while  evaluating  all  of  the  dipole  moment 
matrix  elements  to  account  for  the  electrical  non-linearity.  In  Table  VI,  all  absorption 
bands  with  significant  intensity  are  listed.  Now  v j  is  only  five  times  more  intense  than 
uz,  in  much  better  agreement  with  experiment.  The  increased  intensity  in  u 3  comes 
largely  from  the  non-linear  nCBB  term  in  the  dipole  moment  function. 

When  the  molecule  is  heated  to  3000K,  many  overtones  appear  due  to  the  increased 
sampling  of  the  mechanical  anharmonicity.  The  calculated  absorption  band  shown  in 
Figure  4  is  quite  complicated  but  may  be  assigned  by  performing  calculations  with  some 
components  of  the  dipole  profile  removed.  With  only  the  perpendicularly  polarized 
intensity  present  (only  dfi/dq2  non  zero)  three  bands  result.  These  all  must  have  odd 
quanta  in  the  bend  mode  w2.  At  low  temperatures  the  two  high  frequency  bands  are 
symmetrically  distributed  around  i/j  and  are  clearly  Uj  +i/2  and  ui-i/2.  As  the  tempera¬ 
ture  increases  both  bands  move  to  lower  frequency  and  their  symmetric  location  about 
v j  is  lost,  as  shown  in  Figure  4.  Also  shown  in  Figure  4  is  the  intensity  pattern  when 
only  dfi/dqi  is  non  zero.  Three  bands  again  result,  and  these  must  all  have  even  quanta 
in  the  bend  quantum  number.  They  are  assigned  as  the  CH  stretch  fundamental  and 
the  overtones  2u2  and  6i/2.  As  the  3:1  Fermi  resonance  is  strong  compared  to  the  6:1 
resonance,50  2u2  is  expected  to  be  more  intense  than  6 1/2,  as  observed.  4v2  is  very  weak 
as  it  is  both  poorly  resonant  (4:1)  and  widely  separated  in  energy. 

The  previous  calculations  are  performed  using  several  approximations.  One  way  to 
test  for  their  applicability  is  to  calculate  the  spectra  using  an  independent  technique. 
Currently,  it  is  not  feasible  to  use  time-independent  methods  at  3000K  due  to  the  very 
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the  normal  co-ordinate  potential,  attributable  the  better  fit  near  the  C2v  region. 

C.  Bond  Angle  Distribution  at  High  Temperature 

The  most  striking  features  of  equations  (4.4)-(4.5)  is  the  multiplicative  factor  of 
sinA0  in  equation  (4.5).  Classically,  a  molecule  with  a  linear  equilibrium  geometry  has, 
at  non-zero  temperature,  zero  probability  of  actually  being  linear,  see  Figure  2.  This 
result  can  be  understood  in  terms  that  linear  molecules  have  two  degenerate  bending 
components.  It  is  quite  likely  that  any  one  of  these  components  be  within  d6  of  0  but 
of  vanishing  probability  that  both  components  will  simultaneously  be  within  this  range, 
as  is  needed  for  the  molecule  to  be  truly  linear.  Dynamical  approaches  can  also  be  used 
to  explain  this  result.  Suppose  that  the  molecule  is  bent  and  given  a  random  velocity  in 
both  directions  perpendicular  to  the  molecular  axis.  This  velocity  will  have  a  component 
in  the  plane  of  the  bent  molecule  and  a  component  perpendicular  to  it.  The  perpendicu¬ 
lar  component  produces  angular  momentum  around  the  molecular  axis  which  must  be 
conserved  as  the  molecule  moves.  If  the  molecule  is  to  become  linear  then  the 
corresponding  moment  of  inertia  would  vanish  and  so  the  rotational  kinetic  energy 
would  become  infinite!  Motion  on  the  true  four  dimensional  surface  can  be  viewed  as 
motion  on  a  reduced  three  dimensional  surface  with  a  bend  potential  modified  to  include 
this  centrifugal  barrier  which  forces  the  trajectory  to  stay  out  of  the  linear  region. 

Quantum  mechanics  modifies  this  purely  classical  picture  primarily  by  broadening 
the  distributions  due  to  zero  point  motion.  The  width  of  the  ground  state  wave  function 
is  about  15°  and  so  the  classical  probability  functions  are  only  slightly  broadened. 
These  effects  are  included  in  the  subsequent  spectral  calculations:  indeed,  wave  packet 
dynamics  can  be  used  to  construct  the  quantum  probability  functions,  if  desired.  For 
stiff  molecules  like  HCN  at  reasonable  temperatures,  the  peak  in  this  probability  func¬ 
tion  is  close  to  180  *  and  the  molecule  can  be  treated  as  being  truly  linear.  On  the  other 
hand,  for  floppy  molecules  such  as  C2H  at  high  temperatures,  there  are  large  displace¬ 
ments  from  linearity.  The  classical  average  angle  for  C2H  is  145  ’  at  3000K  and  as  a 
result  it  should  not  be  regarded  as  a  linear  molecule  at  this  temperature.  This  effect 
dominates  the  A(2A')»-  X  spectrum. 

D.  The  Intersection  Region  of  the  Surfaces 

In  the  linear  configuration  the  X  and  A  surfaces  cross  when  Tcc  =  1.348A  at  an 
energy  of  3870cm-1.  Shih,  Peyerimhoff  and  Buenker19  obtained  similar  results.  This 
crossing  is  very  close  to  the  minimum  energy  configuration  of  the  excited  state  and  its 
ramifications  are  expected20, 47  to  be  severe  on  low  temperature  spectra.  When  the 
molecule  bends,  the  crossing  of  the  X  and  A(2A')  surfaces  becomes  sharply  avoided.  This 
is  shown  in  Figure  3,  where  the  vertical  excitation  energies  from  the  X  surface  to  both 
the  A(2A')  and  A(2A")  surfaces  are  plotted  as  contours  at  A rCH  =  0  (  rCH  =  1.056A  ). 
The  value  of  AfCc  at  the  crossing  of  the  X  and  A(2A")  states  is  a  relatively  insensitive 
function  of  bond  angle,  however,  the  crossing  of  the  X  and  A(2A')  states  is  quite  sensi¬ 
tive  to  the  bond  angle.  Because  of  the  rapid  increase  in  energy  of  the  A(2A  )  state  as  the 
molecule  bends,  the  intersection  point  moves  to  very  large  CC  distances.  Note  that  the 
potential  contours  are  somewhat  insensitive  to  the  value  of  Arc„  at  which  they  are 
evaluated. 

Also  plotted  in  Figure  3  are  contours  of  equal  probability  for  the  VJxl(Ar)  poten¬ 
tial  at  3000K.  These  contours  are  centered  around  a  bond  angle  of  145  0  ,  as  noted  previ¬ 
ously.  The  overlap  of  this  set  of  contours  with  the  previous  set  demonstrates  that  the 
intersection  region  is  not  significantly  populated  at  3000K.  Not  only  is  it  improbable  to 
get  direct  Franck-Condon  transitions  into  this  region,  but  also  the  centrifugal  barrier 
prevents  the  excited  state  trajectories  from  entering  it.  Thus,  the  non-adiabatic  effects 
are  reduced20,47  as  the  temperature  increases. 


-  12- 


should  not  be  run  on  these  surfaces  at  energies  near  either  the  saddle  point  energies  or 
end-over-end  rotation  energies.  Figure  2  is  used  to  demonstrate  a  posteriori  that  such 
trajectories  are  not  common  at  a  temperature  of  3000K.  It  shows  the  classical  ensemble 
averaged  values  of  the  bond  lengths  and  angles  at  3000K  for  both  the  Vjx)(Ar)  and 
Vjxl(q)  potentials,  and  these  show  little  penetration  into  the  poorly  represented  regions. 
Expressions  for  these  average  values  can  be  obtained  by  considering  the  partition  func¬ 
tion,  which  for  the  normal  co-ordinate  potential  is 

Qn  =  f  d<t\dq2,dqudqz  e  ‘(q)  .  (4.1) 

«/  —  00 

Noting  that  Vjsi(q)  are  only  functions  of  q2  =  (92«2  +  ?2fc2)  »  we  can  introduce  the 
"polar”  co-ordinates  q2  and  <t>,  as  is  done  in  equation  (3.2),  simplifying  the  partition  func¬ 
tion  to 

Qn  =  2*  d1sf0  dl2  92  «  ^  ■  (4-2) 

A  similar  result  applies  to  the  partition  function  expression  written  in  terms  of  valence 
coordinates,  resulting  in 

Qv  =  foB°drcHdrccfo*d0  sinA*  e^V'IX,(Ar)  .  (4.3) 


Here,  the  sin  AS  term  arises  from  the  transformation  from  cartesian  to  valence  co¬ 
ordinates,  equation  (3.1).  The  probabilities  of  being  within  drCH  of  TCH ,  dT cc  of  Tcc 
and  d  0  of  9  are  given  by 
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respectively.  Parallel  equations  exist  for  these  probabilities  written  in  terms  of  normal 
coordinates. 

In  Figure  2  the  bond  length  distributions  are  well  within  the  reliability  range  of  the 
fitted  X  potential  functions.  Unfortunately  the  angle  distribution  has  a  slight  tail  to 
large  angles,  into  the  end-over-end  hydrogen  rotation  region.  Trajectories  which  start  in 
this  region  of  the  potential  are  not  run  in  the  subsequent  spectral  calculations.  This  is 
not  a  significant  error  in  the  calculation.  It  is  clear  that  these  potential  surfaces  are 
inappropriate  at  temperatures  significantly  higher  than  3000K  for  which  end-over-end 
rotation  needs  to  be  included  explicitly. 

A  method  for  estimating  the  reliability  of  the  A  surfaces  is  to  take  the  most  prob¬ 
able  regions  on  the  X  state  and  look  at  the  excited  state  surfaces,  not  only  in  this  region, 
but  in  all  regions  energetically  accessible  to  excited  state  trajectories  starting  with  these 
ground  state  geometries.  The  excited  state  data  points  are  chosen  to  span  this  entire 
region,  and  about  1%  of  all  the  trajectories  which  were  run  were  found  in  poorly 
represented  regions. 

The  probability  distribution  functions  for  both  the  normal  and  valence  co-ordinate 
potentials  are  given  in  Figure  2  as  a  final  check  on  the  accuracy  of  the  functions. 
Clearly,  both  potentials  are  very  similar.  The  major  differences  in  the  probability  func¬ 
tions  are  first,  a  shift  of  0.015A  in  the  TCH  probability,  attributable  to  inadequacies  in 
the  bend-CH  stretch  potential  (see  Figure  1);  and  second,  the  longer  tail  to  low  angle  on 


Figure  1.  Potential  Contours  in  1000cm-1  above  each  states  own  minimum. 
Top  panel  gives  normal  co-ordinate  contours  at  Bottom  panel  gives  valence 

co-ordinate  contours  at  &TCh  =0. 
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The  calculated  X  state  CH  stretch  harmonic  vibrational  frequencies  are  nearly  the 
same,  3645cm-1  and  3623cm-1,  for  Vjxl/Al*)  and  l;JXl(q),  respectively.  When  anhar- 
monicity  is  included  the  frequency  for  V„'x!(q)  decreases  to  3484cm-1,  slightly  less  than 
the  value  of  3612cm-1  observed13  by  Jacox.  In  the  A  state  the  CH  stretching  frequency 
is  slightly  reduced.  The  CC  harmonic  stretch  frequency  is  calculated  to  be  near 
1975cm-1  in  the  ground  state.  Anharmonic  corrections  reduce  this  to  1900cm-1,  not  far 
from  the  observed18  value  of  1848cm-1.  In  the  A  state  the  corresponding  harmonic  fre¬ 
quency  is  1655cm-1,  and  anharmonic  corrections  will  likely  bring  this  very  close  to  the 
observed12  value  of  1560i  40cm-1.  For  the  bending  vibration  the  calculated  harmonic 
frequencies  are  around  494cm-1:  anharmonic  corrections  reduce  this  to  470cm-1.  This 
frequency  has  never  been  directly  observed.  Carrick,  Merer  and  Curl15  deduced,  by 
analyzing  the  /-type  doubling  in  their  electronic  spectra,  that  it  lies  in  the  range  390- 
475cm-1:  arguing  that  this  range  is  conservative,  they  then  went  on  to  postulate  that 
the  frequency  is  actually  375cm-1.  Clearly,  further  experimental  studies  need  to  be  per¬ 
formed  in  order  to  measure  this  frequency  more  reliably.  In  the  excited  A(2A")  surface 
this  frequency  decreases  slightly  but  it  more  than  doubles  in  the  excited  A(2A')  surface. 
This  is  brought  about  by  the  repulsion  between  the  lower  X  (2A')  and  upper  A(2A')  sur¬ 
faces,  and  has  a  very  important  effect  on  high  temperature  spectra. 

B.  Useful  Range  of  the  Potentials 

Figure  1  contains  contour  plots  of  sections  cut  through  the  surfaces  and  indicates 
the  magnitude  of  the  anharmonicity  of  C2H.  The  top  three  inserts  in  this  figure  display 
Vjsi(q)  cut  through  the  section  q 3  =  0,  while  the  lower  inserts  display  VjSl(Ar)  cut 
through  the  section  A Tcc  =  0.  Here  the  contours,  in  1000cm-1,  are  given  with  reference 
to  the  minima  of  the  respective  potential  surfaces  S.  As  the  A(2IT)  surface  is  displaced 
in  the  CC  stretch  variable  the  lowest  energy  in  both  of  the  A  sections  is  2027cm-1  above 
the  A  state  potential  minimum.  From  the  shape  of  the  contours  we  see  both  the  intrinsic 
advantage  of  writing  the  potential  in  valence  co-ordinates  and  the  effectiveness  of  the 
non-linear  transformation,  equations  (3.20)-(3.22). 

These  plots  extend  in  geometry  far  beyond  the  region  of  the  potential  that  is  fitted 
to  the  ab  initio  data  /most  of  the  data  points  have  energies  <  10000cm-1,  very  few  have 
energies  >  16000cm-  ).  In  the  fitted  regions  the  potentials  fit  the  data  well.  Spurious 
saddle  points  exist  outside  this  region,  however.  Consequently  these  surfaces  are  not 
suitable  for  extrapolation  purposes.  Included  in  Table  V  are  the  geometries  and  energies 
of  the  lowest  found  saddle  point  for  each  potential  surface. 

The  valence  co-ordinate  potentials  are  slightly  less  susceptible  to  low  lying  saddle 
points.  In  some  respects  this  is  a  disadvantage  as  the  normal  co-ordinate  potentials 
represent  better  the  turn-over  of  the  bend  potential  which  occurs  as  the  C2v  geometry  is 
approached.  From  Tabte  V  the  valence  potential  predicts  that  a  saddle  point  lies  on  the 
X  surface  at  Tc„  =  1.33A  ,  Tcc  =  1.34A  ,  0  =  57  *  and  an  energy  of  14789cm-1.  As  the 
largest  angle  used  in  the  CASSCF  data  is  100  ’ ,  this  point  represents  a  considerable 
extrapolation.  Subsequent  to  the  fitting  of  the  potentials,  more  CASSCF  calculations 
have  been  performed  and  the  C2v  transition  state  for  end-over-end  hydrogen  rotation  cal¬ 
culated  to  be  at  T CH  =  1.28A  ,  T cc  ~  1.29A  ,  0  =  63  '  and  an  energy  of  14509cm-1.  As 
the  CASSCF  energy  increases  quite  slowly  when  all  bonds  are  simultaneously  extended, 
the  accuracy  of  the  extrapolation  of  the  norma)  co-ordinate  potential  is  quite  remarkable. 
On  the  other  hand  the  valence  co-ordinate  potentials  tend  to  climb  rapidly  at  small  0, 
showing  no  sign  of  turning  over.  If  greater  accuracy  is  needed,  then  a  potential  written 
in  terms  of  the  sine  of  the  angle  between  the  CC  bond  vector  and  the  vector  joining  the 
hydrogen  atom  to  the  mid-point  of  the  CC  bond  vector  should  be  used.  Such  a  function 
would  possess  the  correct  angle  boundary  conditions. 

At  moderate  to  high  energy  the  surfaces  are  quite  ergodic:  a  trajectory  will  quickly 
sample  most  of  phase  space  that  is  energetically  accessible  to  it.  Thus  trajectories 
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TABLE  VI.  Vibrational  Absorption  Frequencies  and  Intensities  from  the  X 
state  Normal  Co-ordinate  Porential. 


Mode 

Component 

Frequency 

(cm-1) 

Relative 

Intensity 

v2 

n 

1. 

2 1/2 

E4 

980 

^3 

E4 

0.17 

n 

2326 

I 

2^2  +  uz 

E4 

2857 

0.03 

E4 

3484 

0.94 

2 1/3 

E4 

3778 

n 

3940 

^1+^2 

E4 

5380 

2i/j 

E4 

6855 

L 


12 


57478 


ordinates.  The  functional  forms  used  to  represent  these  components  are 

=  Mo  +  Ph  9i  +  Me  93  +  Pm  9  ]2  +  Phc9\9z  +  Pcc  9 z  +  Pbb  92  +  (3-20) 

2  2  22  2  22  A 

Mwi  9 1 9 2  +  Pcbb  9s9i  +  Phhbb  9 i  92  +  Phcbb  9 1 9 39 2  "+  Pccbb  9 3  ?2  +  Pbbbb  9 2 

and 

Mo  =  ?2  (  Pb  +  Pub  9i  +  Pcb  9z  +  Phhb  9*  +  Phcb  9\9i  +  Pccb  9  3*  +  Pbbb  92  +  (3.21) 

S  2  2  3  2  2  \ 

Phhhb  9 1  +  Phhcb  9]  ?3  +  Phccb  9] 9 3  +  Pcccb  9 3  +  Phbbb  9i92  "1"  Pcbbb  9  s? 2  )  • 


The  electronic  transition  A(2II)<-  X(2S+)  has  different  symmetry  properties.  It  is  per¬ 
pendicularly  polarized:  fit  is  used  to  represent  these  components.  When  the  molecule 
bends,  the  in-plane  component  of  the  transition  moment  is  associated  with  the 
A(2A')<-  X  transition  and  the  out-of-plane  component  is  associated  with  the  A(2A")<-X 
transition.  Bend  displacements  also  introduce  some  parallel-polarized  vibronic  intensity: 
fi0  is  used  to  represent  this  component  of  the  transition  moment.  This  intensity  is  asso¬ 
ciated  with  the  A(2A')«-  X  transition  only. 


Minimization  of  the  root  mean  square  error 

Vz 

Erm,  =  —  £  (Passer -pm  )2 
"  » 


(3.22) 


is  the  condition  used  to  determine  the  parameters  in  these  expansions.  As  the  number 
of  ab  initio  data  points  for  the  transition  moments  is  quite  small,  some  of  the  cross 
terms  in  these  expansions  are  constrained  to  be  zero.  Different  refinements  are  per¬ 
formed  for  both  of  the  odd  and  even  functions  and  for  both  the  pure-stretch  and  bend- 
containing  subsets  of  the  functions  The  calculated  constants  are  given  in  Table  III,  and  a 
summary  of  the  fits  is  given  in  Table  IV  including  the  error  term  Ermi ,  the  number  of 
free  parameters,  and  the  number  of  data  points.  More  details  are  available  from  the 
authors.  Note  that  the  bending  X  dipole  moment  functions  fit  the  data  points  with 
bond  angles  ■:  145’  very  well:  most  of  the  errors  reported  in  Table  IV  arise  from  the 
data  points  with  bond  angle  in  the  range  100 '  to  145  * . 

These  expansions  in  terms  of  the  dimensionless  normal  co-ordinates  converge 
quickly.  The  only  high-erder  constants  that  are  appreciable  are  PbbiPcbb  and  Pcb • 
Similar  fits  to  the  «6  tntl  '  data  have  been  attempted  by  expanding  in  powers  of  the 
valence  co-ordinates  AT  i ... ,  »re  not  reported  here  as  the  fits  are  unstable,  having 
large  high  order  terms  Apparently  the  normal  co-ordinates  provide  a  better  description 
of  the  dipole  moment  surface  than  do  the  valence  co-ordinates. 

IV.  PROPERTIES  OF  THE  POTENTIALS 
A.  Vibration  Frequencies 

Table  V  gives  details  of  a  harmonic  analysis  of  the  potential  minima.  All  minima 
have  linear  geometries  and  their  bond  lengths  may  be  obtained  directly  from  the  param¬ 
eters  given  in  Table  1.  To  get  an  indication  of  the  effects  of  anharmonicity,  the  vibra¬ 
tional  hamiltonian  matrix  for  the  Vjxl(q)  potential  is  constructed  in  the  norma)  co¬ 
ordinate  basis,  truncated  at  states  with  more  than  a  total  of  eight  quanta  (a  grand  total 
of  495  states  are  included),  and  then  diagonalized.  A  product  form  Gauss-Hermite  qua¬ 
drature  is  used  to  evaluate  the  matrix  elements.  Some  of  the  eigenvalues  are  given  in 
Table  VI,  relative  to  the  zero  point  energy  which  is  found  to  be  3240cm'5.  The  one  and 
two  quanta  levels  are  believed  to  be  converged  to  ±  10cm'5  and  ±  40cm"5,  respectively: 
note,  however,  that  no  rotational  corrections^7  are  included  in  these  calculations. 


3000K 


1000K 


Each  line  is  broadened  by  unresolved  hot  bands  ]2q  3  expected  slightly  to  the  red  of  the 
main  ^  3  lines.  There  is  also  a  small  intrinsic  broadening  due  to  the  use  of  the  frozen 
gaussian  approximation  for  propagation  on  an  anharrnonic  potential  surface.26  At  3000K 
the  spectrum  contains  the  cold  lines  q  3  as  well  as  two  hot  band  progressions  "3  and  \  3. 
Lines  are  drawn  on  Figure  5  to  highlight  these  hot  band  progressions.  The  3000K  spec¬ 
trum  is  broadened  due  to  the  large  number  of  lines  present,  and  only  partial  structure 
remains  at  60cm-1  resolution. 

The  calculated  spectra  for  the  A(2A")<-  X  transition,  Figure  7,  are  much  more  diffi¬ 
cult  to  interpret.  There  is  both  a  perpendicular  band  based  upon  the  allowed  origin  O 
and  also  a  parallel  band  based  upon  the  vibronic  origin  j2.  The  perpendicular  band  at 
OK  is  almost  identical  to  the  perpendicular  band  for  the  A(2A")*-  X  transition  except 
that  the  zero  point  correction  shifts  the  origin  up  to  3960cm-1.  At  1000K  the  perpendic¬ 
ular  spectra  of  the  two  states  are  quite  different.  Hot  bands  nn  2  are  spaced  progressively 
at  550cm-1  to  the  blue  of  all  of  the  OK  lines  and  fill  the  gaps  between  the  q  3  progression. 
The  intensities  of  these  lines  are  highly  non-Boltzmann  with  both  \2  and  |2  being 
stronger  than  This  enhancement  comes  from  the  high  anharmonicity  of  the  excited 
state  bend.  The  A(2A')  potential  decreases  markedly  in  curvature  at  large  angle,  thus 
increasing  the  overlap  between  its  wave  functions  and  the  ground  state  wave  functions. 

There  are  problems  with  the  frozen  gaussian  approximation  which  even  at  OK  have 
more  profound  effects  than  just  missing  the  changes  in  the  zero  point  levels.  Even 
though  there  is  no  displacement  in  the  bend  co-ordinate  between  the  minimum  energy 
configuration  of  the  ground  and  excited  states,  Frank-Condon  factors  linking  levels  like 
o2  can  be  significant.  This  is  due  to  the  large  frequency  change  between  the  ground  and 
excited  state  bending  vibrations.  Thus  a  second  progression,  based  upon  the  origin  O 
and  spaced  at  2100cm-1,  should  be  present.  Thawed  gaussians  know  about  the  fre¬ 
quency  change  and  trial  calculations  performed  using  them  indicate  that  is  about  five 
times  weaker  than  ^2.  Thus  the  error  in  using  frozen  gaussians  is  not  great.  At  higher 
temperatures  there  is  motion  in  the  bend  co-ordinate  and  thus  the  calculation  knows 
about  the  frequency  change.  In  the  1000K  spectrum  the  shoulder  at  5960cm-1  to  the 
left  of  the  fifth  strong  band  is  ^2  and  a  progression  in  i/3  can  be  seen  based  upon  this 
line.  The  large  number  of  lines  resulting  at  high  energy  can  not  be  distinguished  at 
60cm-1  resolution. 

Results  for  the  forbidden  parallel  band,  also  given  in  Figure  7,  show  essentially  the 
same  features.  The  problem  with  the  frozen  gaussian  approximation  at  OK  in  not 
detecting  the  true  excited  state  bend  frequency  now  results  in  the  vibronic  origin  $ 
being  seen  550cm-1  to  the  red  of  where  it  should  be.  Again  at  1000K  the  lines  move  to 
their  correct  location. 

Figure  8  gives  the  total  A«-  X  spectra  at  OK,  1000K  and  3000K  as  a  sum  of  the 
individual  A(2A')<-X  and  A(2A")»-  X  spectra.  The  doublets  in  the  OK  spectra  arise  due 
to  the  difference  in  the  zero  point  energies  of  the  two  A  state  components.  Inadequacies 
in  this  spectrum  have  been  discussed  previously:  the  red  shift  of  550cm-1  in  the  parallel 
transitions  and  the  absence  of  the  weak  3  2  progression.  At  1000K  the  large  number  of 
new  lines  which  appear  washes  out  the  structure  at  high  energy.  This  continues  as  the 
temperature  rises  and  at  3000K  very  little  structure  remains  at  60cm-1  resolution. 
There  is  a  large  blue  shift  in  the  intensity  distribution  at  high  temperatures  due  to  the 
repulsive  nature  of  the  A(2A')  surface  combined  with  the  large  value  of  the  average 
ground  state  bend  angle.  The  maximum  intensity  in  the  spectrum  is  shifted  to 
9500cm-1.  Some  indication  of  a  two  band  structure  is  present  in  the  spectra  with  the 
A(2A")«-X  band  appearing  as  a  shoulder  4000cm-1  to  the  red  of  the  intense  A(2A')<-X 
band. 


3000K 


1000K 


VI.  CONCLUSIONS 

This  work  is  the  first  application  of  the  high  temperature  spectral  formalism  of  Rei- 
mers,  Wilson  and  Heller26  to  either  infrared  spectra  or  the  spectra  of  polyatomic 
molecules.  The  method  is  seen  to  be  both  simple  and  practical,  requiring  little  more 
effort  than  pure  classical  mechanics.  Some  problems  appear  because  of  the  use  of  frozen 
gaussians  and  classical  temperatures  but  these  fade  in  significance  as  the  temperature  is 
increased  and  the  density  of  spectral  lines  increases.  Also  at  high  temperatures  the  dom- 
inately  bent  character  of  the  molecule  reduces  the  effects  of  the  Born  Oppenheimer 
approximation,  allowing  the  A<-  X  absorption  to  be  represented  as  the  sum  of  the 
A(  A')<-X  and  A(2A")«-  X  absorptions. 

The  C2H  potential  surfaces  have  distinctive  properties  which  are  not  characteristic 
of  typical  laboratory  molecules  like  HCN.  First  C2H  is  floppy,  leading  to  large  thermal 
bend  displacements.  When  this  couples  with  the  appreciable  mechanical  anharmonicity 
and  electrical  non-linearity,  many  overtone  and  combination  bands  appear  in  the 
infrared  spectra.  Both  the  natures  of  the  Fermi  resonances  involved  and  the  energy  level 
differences  are  seen  to  be  important  in  redistributing  intensity  from  the  fundamentals  to 
the  overtones.  In  particular  the  CC  stretch  fundamental  is  predicted  to  get  most  of  its 
intensity  from  anharmonic  effects,  and  its  intensity  is  predicted  to  be  highly  temperature 
dependent.  The  A(2A")<- X  electronic  transition  is  relatively  simple  but  the  A(2A')«-  X 
absorption  is  complicated  by  the  large  thermal  motion  and  bend  frequency  disparity. 
We  see  many  blue  shifted  hot  bands  with  highly  non-Boltzmann  intensity  distributions, 
and  progressions  in  a  mode  that  is  not  displaced.  The  center  of  the  absorption  band 
shifts  4000cm'1  to  the  blue  when  C2H  is  heated  from  OK  to  3000K.  Its  intensity 
increases  considerably  due  to  the  forbidden  nature  of  the  parallel  transition,  leaving  the 
unaffected  A(2A")*-  X  band  as  a  shoulder  to  the  red  of  the  A(2A')«-  X  band.  A  very 
large  number  of  lines  appear  in  the  hot  spectra  and  these  blend  into  a  continuum  at 
60cm'1  resolution. 

Experimentally  only  the  region  3500-4200cm-1,  in  the  tail  of  the  absorption  band, 
has  been  investigated.15  Within  this  region  the  calculated  intensity  due  to  the  electronic 
transitions  (Figures  6-8)  is  much  greater  than  the  vibrational  intensity  (Figures  4-5). 
This  supports  the  conclusions  of  Carrick,  Merer  and  Curl15  that  none  of  their  observed 
bands  in  this  region  is  primarily  of  vibrational  nature.  This  work  shows  that  the  overall 
band  structure  contains  a  lot  of  information  about  the  potential  surfaces  and  therefore  it 
is  desirable  that  the  entire  band  be  experimentally  observed  in  the  gas  phase. 


We  have  shown  that  it  is  possible  to  calculate  ab  initio  the  often  complicated 
spectra  of  small  molecular  fragments  like  C2H.  There  exists  many  chemical  species  in 
this  category,  ranging  from  covalent  second  row  type  molecules  to  largely  ionic  molecules 
like  NaOH  and  NaCl  to  fragments  of  transition  metal  complexes  and  clusters.  In  general 
these  are  difficult  species  to  observe  and  characterize  in  the  gas  phase  so  that  their  spec¬ 
tra  are  difficult  to  obtain.  Molecules  of  this  type  are  important  though  as  reaction  inter¬ 
mediates  and  in  astrochemistry,  and  they  are  often  found  at  high  temperature.  It  is 
hoped  that  the  CASSCF  and  wave  packet  techniques  used  here  will  be  of  general  use  for 
predicting  their  spectra. 
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